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ABSTRACT 



We discuss a new solution of the Kompaneets-equation for physical situations in which low frequency photons, forming relatively 
narrow spectral details, are Compton scattered in an isotropic, infinite medium with an intense ambient blackbody field that is very 
close to full thermodynamic equilibrium with the free electrons. In this situation the background-induced stimulated Compton scatter- 
ing slows down the motion of photons toward higher frequencies by a factor of 3 in comparison with the solution that only takes into 
account Doppler-broadening and boosting. This new solution is important for detailed computations of cosmic microwave background 
spectral distortions arising due to uncompensated atomic transitions of hydrogen and helium in the early Universe. In addition we 
derive another analytic solution that only includes the background-induced stimulated Compton scattering and is valid for power-law 
ambient radiation fields. This solution might have interesting applications for radio lines arising inside of bright extra-galactic radio 
source, where according to our estimates line shifts because of background-induced stimulated scattering could be amplified and even 
exceed the line broadening due to the Doppler-effect. 
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1. Introduction 

The recombination of He III (4500 < z < 7000) in the 
Universe produces specific quasi-periodic distortions of cos- 
mic microwave background (CMB) in the Rayleigh-Jeans 
part of the CMB spectrum , at dm and cm wavel ength 
dDubrovich & Stolvarovl Il997t iRubino-Martin et all |2007|) . It 
is also well-known that any energy release below z ~ 
5 x 10 4 will lead to a characte ristic broad y-type dis- 
tortion dZeldovich & Sunvaevl [T969) of the background ra- 
diation. In this case bec ause of uncompensated loops 
dLvubarskv & Sunvaevlll983l) . trying to restore full equilibrium, 
similar features should originate from even higher redshifts 
(Chluba & Sunvaev 2008). However, their amplitude and phase- 
dependence should differ from those of the distortions arising 
during He III recombination. 

It is clear that in both cases the multiple scattering of pho- 
tons by hot electrons should lead to broadening and shifting 
of these spectral features. Here we investigate the evolution of 
narrow lines due to both stimulated scattering in the presence 
of the much brighter CMB background with temperature T y = 
2.7 (1 + z) K, and Doppler-broadening and boosting. According 
to our analytical solutions, in this situation the background- 
induced stimulated Compton scattering at low frequencies slows 
down the motion of photons toward higher frequencies by a fac- 
tor of 3 in comparison with the solution t hat only takes into ac- 
count Doppler-broadening and boosting dZeldovich & Sunvaevl 
[l969h . 

We present two analytical solutions demonstrating this fact. 
In the first we only include the pure effect of stimulated scatter- 
ing in the presence of a bright power-law ambient photon field, 



neglecting the Doppler-effect (see Sect. [3), while in the second 
solution we take into account both effects simultaneously, but 
restrict ourselves to the case of a low-frequency blackbody am- 
bient photon field with temperature T y ~ T e , where T e denotes 
the electron temperature (see Sect. HJ. The second solution also 
clearly shows that the broadening of the weak lines in this situa- 
tion depends only on the y-parameter defined by T e , even though 
the evolution of the ambient CMB blackbody spectrum itself is 
described by y oc r e - T y . 



2. Summary of previous analytic solutions to the 
Kompaneets-Equation 

The repeated Compton scattering of photons by thermal elec- 
trons in isotropic, infinite med ia can be described using the well- 
known Kompaneets-equation dKompaneets| [l956) 



dn(x e ,y e ) _ 1 d 
dy e x\ dx e 



dn 

dx. 



+ n(l + n) 



(1) 
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where n(x e ,y e ) = c 2 I v /2hv 3 is the photon occupation number, 
y e = f -^r N e cr T cdt is the Compton y-parameter, where T e de- 
notes the electron temperature, and x e - hv/lcT e is the dimen- 
sionless frequency. 

As is well understood, the first term in brackets (oc d^n) de- 
scribes the diffusion of photons along the frequency axis, the 
second term (oc «) the motion of photons towards low frequen- 
cies due to the reco;7-effect, and the last term (oc « 2 ) the effect of 
stimulat ed scattering, which physi cally is also related to recoil 
(e.g. see lSazonov & Sunyaevll2000l) Eq. (ffli has been studi ed in 
great detail, both numerically (e.g. iPozdniakov et al.ll 19831) and 
analytically in several limiting cases, but, to our knowledge, no 
general analytic solution was found. 
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In situations when the recoil-term and stimulated scatterings 
are no t important (i.e. n[n + 1] <K d r _ri). IZeldovich & Sunvaevl 
(1969) gave a solution for arbitrary value of y e , which reads 



n(v,y e )\ m = 



^ftny^ 



I 



_ (ln[v/v]+3y c ) 2 

n (v)e 4 * 



dv 



(2) 



Here hq(v) is the photon occupation number at frequency v and 
y e = 0. For y e <K 1 the broadening of an initially narrow line 
is given by Av/v ~ ±2 Jy~ e In 2. Due to Doppler-boosting, the 
maximum of the specific intensity I v moves along vn ne (y e ) = 
voe 3>e , so that for y e <k 1 one has Av/v ~ 3y e , implying that 
photons are upscattered. 

In the case when the d iffusion term and stimulated scat- 
terings are not important. | Aronsl (Il97lh and independently 
llllarionov & Svunvaevl ( 1 19721) . using different mathematical ap- 
proach, gave a solution which can be written in the form 



"0,T T )| r ecoil = 



1 



. no i | . 

[1 — COTj] 4 \1— 0>Tt/ 



(3) 



where u> = and t t 



N e crjct is the Thomson-scattering 
optical depth. This solution simply describes the motion of pho- 
tons towards lower frequencies, where for initial photon distri- 
bution, «o(v) = A6(v - vo)/v 2 , at later time the line is located 
at viine(Tx) = vo/[l + woTj]. For wtt «C 1 the line shift due to 
recoil is given by Av/v ~ -cotj. This shows that at low frequen- 
cies (hv <§: kT e ) the recoil-effect can be neglected in comparison 
with Doppler-boosting. 

The third analytic solution was found i n the case when only 
the st i mulated scatte ring term is important dZeldovich & Levichl 
1968; Sunvaev1 ll970h . Here the solution is determined by the 
implicit equation 



0(/)-2- 



(4) 



with /(v, tt) = v 2 n(v, tt) and where <f>(z) can be found from 
the initial condition (</>(z) = f l (z), where f l (z) is the inverse 
function of f(v, Tt) at Tt = 0). 

Note that the last two solutions do not depend on the temper- 
ature of the electrons, since in both cases physically the motion 
of the electrons in lowest order does not play any role. 

3. Solution for purely background-induced 
stimulated scattering 

In this Section we discuss the situation when photons from a 
low-frequency spectral feature are scattering off thermal elec- 
trons at temperature T e within an intense photon background that 
has a power-law spectral dependence, i.e. «back(i / ) = Aback v$ » 
1 . Here «o,back = Ab ac kvjj » 1 is the occupation number of back- 
ground photon field in the vicinity of the spectral feature at vo, 
and p denotes the power-law inde?|3. 

If we neglect the evolution of the background photon fielcQ 
and assume that recoil and Doppler broadening are negligible, 
then, inserting n = «back + An into Eq. ((1), one is left with 



<9An(v,r T ) 2 h d 



<9r T 



-7 7 tt v «backO) An . 



V OTpC 



2 dv 



1 When the spectral intensity scales like /„ oc v a , one has /? = -a— 3. 
For the Rayleigh- Jeans part of a blackbody spectrum /3 = - 1 . 

2 This assumption is for example justified when the frequency shifts 
of the background radiation field are relatively small. 



For initial photon distribution Ano(v) = An(v, tt = 0) the solu- 
tion is given by 



Anr 



An(v,T T )lback-ind = 



[1 -2(]3+ l)n b ack(v)wT T p 



(6) 



This solution describes the purely background-induced motion 
of photons along the frequency axis. Due to the large factor 
2 nback s> 1 the speed of this motion is increased, so that even 
for ujtt <k 1 the line shift can still be considerable. For fi > — 1 
this solution is valid only at v < l/[2(/3 + 1) A bac k ^t] 1W+1) - 
For initial photon distribution, Ano(v) = A 6(v - vo)/v 2 , the 

line is located at v*(tt) = vo/[l+2(j8+l) nb ac k(vo) woTt]^- Here 
the solution works for all /J > - 1 , while for fi < - 1 it is valid 
only for v > [2\f3 + 1| A back For the condition 2(J3 + 

l)«backOo) woTt <k 1 one has Av/v ~ -2 n bac k(vo) wo t t . This 
shows that independent of /3 photons are always moving toward 
lower frequencies, where the speed of this motion is increased 
by 2 nback » 1 as compared to the recoil-case (see Sect.|2]). It is 
also remarkable that for ft = -1, i.e. nb ac k 1/v, one has Av/v * 
-2 Aback ^72 , so that the line shift is independent of frequency. 

It is also possible to rewrite this expression in terms of the 
brightness temperature 



2 rback 



?b(v) = 



hv 



2kv 2 



«back(v), 



(7) 



of the background field in the vicinity of the spectral feature, 
yielding 



Av 
v 



back-ind 



„ kT b (v) 

-2 5-TT. 



(8) 



For 7b > jT e , it is clear that the line shift due to background- 
induced stimulated scattering exceeds the shift towards higher 
frequencies due to Doppler-boosting. Therefore the net motion 
of the line center can be directed towards lower frequencies, even 
if one includes the Doppler-term. Note that for a Rayleigh- Jeans 
spectrum r b is equal to the physical temperature of the radiation 
field. 

Comparing with the case of pure Doppler-broadening (see 
Sect.©, it is also clear that fovy e <k 1 andy e > 4 [T e /T b (y)] 2 In 2 
the background-induced line shift becomes larger than the 
FWHM Doppler-broadening connected with the diffusion term. 
In this case one can neglect the effect of line-broadening due 
to the Doppler-term, and, according to Eq. ((H), even for small 
Thomson optical depth can still obtain a rather significant line- 
shift. In radio spectroscopy is possible to determine tiny shifts 
or broadening in the frequency of narrow spectral feature (for 
example of 21 cm lines from high redshift radio galaxies). 



3.1. Rayleigh-Jeans limit (fi--\) 

For/? = -1 Eq. © becomes <9 TT An(v, r T ) = £ d v v 3 An, and 
with the replacement s = v 3 An and £ = ln(v), can be cast in 



(5) the form d Tj s(£, tt) : 



2hA b ., 



■ d^s. Using the characteristic of this 
equation, or by directly taking the limit fi — > — 1 from Eq. ((6), 
one readily finds 



A"(v,T T )tck-ind =«' 



6,1,,, 



"e? An lve ^ , 



(9) 
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so that for Ano(vo) = A 6(v - vq)/v , the line will be centered at 

-2A At t 

v*(tt) = vo e b " ck . If we assume that «back = kTy/hv, i.e. we 
are in the Rayleigh-Jeans part of a blackbody photon field with 

T y — 7b = T e , then we have v*(tt) — V(,e »«< 2 T . This shows 
that due to stimulated scattering of photons within a blackbody 
ambient radiation field are moving with a speed oc 2y e towards 
lower frequencies. This is 3/2 times smaller than the line shift 
due to Doppler-boosting (see Sect. but is directed in the op- 
posite direction. Below we will discuss this case in more detail, 
also including the broadening of the line due to the Doppler- 
effect. 



3.2. Evolution of the background-field and the relative 
velocity of a spectral feature for small line-shifts 

In the derivation given above we have neglected the evolution 
of the background radiation field. This is possible when the line- 
shifts Av/v are small. If we insert n = «b ac k into the Kompaneets- 
equation (Q]i, and assume thaty e is sufficiently small, then we can 
directly write the solution for the change of the photon occupa- 
tion number as 



An ba 



"back 



y e {/3(3 + (3) + (4 + /3)x e + 2«back(2 + fi)x e 



(10) 



With this one can now find the corresponding overall shift of the 
background radiation field at initial frequency x e Solving the 



equation n ba ck(x e ,o) = «back(x e )[l + i 



Mback(-rc) 

for 6 = -Av/v under the condition ft + one finds: 



],withx e = x e fl(l+6), 



Av 
v 



back 



2(2+0)kT b (v) 4 hv 

"Tt + - J T T 



+ (3 + fi + X e ) • 



P m e c 2 
kT e 

— T T T ■ 



(11) 



With |Av/v| <K 1 one can give the condition under which it is 
possible to neglect the evolution of the background field. 

It is important to mention that Eq. ( fTOb does only conserve a 
Wien spectrum (« = e~ Xc ~ 1 for x e <k 1) and Rayleigh-Jeans 
spectrum (n — 1 /x e ) to order x e . This is because higher order 
terms cannot be taken into account by simple power-laws. Only 
n — 1 /[e Xc+ft - 1] with constant /i and T e = T 7 is truly invariant. 

4. Scattering of low frequency photons by free 
electrons in an intense ambient blackbody field 

If we assume that at y e — the initial photon field is given by 
n - n p i(x e )+An, where n v \(x) - 1 /[e x - 1] is the blackbody occu- 
pation number, x e = hv/kT e , and An is a small (An/n p i(x e ) <k 1) 
spectral distortion, then with d Xc n p i(x e ) = -n p i(x e )[l + n p \(x e )] 
and dy c n p i(x e ) = 0, from Eq. ([1) one has 



dAn(x e ,y e ) 1 
dy e x\ 



d 



dAn 

+ An(l + 2« p i(x e ) + An) 

ox. 



(12) 



If we now neglect terms of 0(An 2 ) and assume that x e <k 1, and 
hence n p \(x e ) * 1 /x e » 1, then Eq. ( fT2b reads 



dAn(x e ,y e ) 1 d 



dy e 



xl dx. 



dAn + 2^ n 

(3»^Cp Xq 



(13) 



Introducing the variables ^ = ln(x e ) and s = xl An we then have 
dy c s(^,y e ) = d\s - d^s. Similar to Zeldovic h & Sunvaevl fl969) 



we now transform to z — £ - y e so that this equation reduces to 
the normal diffusion equation d ye s = d\s. Therefore the solution 
ofEq. ©is 



i r v 3 

An(v,y e ) = -rAn(v, 



_ (ln[v/i>]-y c ) 2 j- 

0)e 4 >-< —. (14) 
v 



Not including the term 2^ in Eq. ( fT3l ). simply leads to the re- 
placement y e — * 3y e in the exponential term of Eq. (fT4T >. and af- 
ter absorbi ng the factor v 3 /v 3 - exp(-3 ln[v/vl) one obtains the 
solution of Ze ldovich & Sunvaevl d 1969b in the form of Eq. (f2j). 

If we assume An(v, 0) = A 5{v - vo)/v 2 for the initial distor- 
tion then from Eq. ( TT4b we obtain 



A«(v,y e ) = 



x e 



(ln[v/v ]-v L .) 2 

4>e 



(15) 



Without the inclusion of stimulated scattering in the ambient 
blackbod y field (replacement y e — > 3y e in Eq. (fT4b). as in the 
case of IZeldovich & Sunvaevl (1 19691) . one finds 



An zs (v,y e ) = 



V47!- .Ve V 3 

A 



(ln[v/v ]-3v e ) 2 

4, C (16a) 

(ln[v/v ]+3v c ) 2 



(16b) 



Comparing Eq. dT3T > and ( I16ab one can see that in terms of energy 
(A7 V oc v 3 An) the maximum of the distribution is moving like 
VmaxCye) = vo e y °, when including the effect of stimulated scatter- 
ing, while it moves like v^ x (y e ) = vo e 3v » without this term. On 
the other hand, in terms of photon number (AA^ oc v 2 A«), from 
Eq. (TT3T > and ( I16al ) one finds that the maximum of the distribu- 
tion is moving like v max (y e ) = Voe~ y °, when including the effect 
of stimulated scattering, while it moves like v^j x (j e ) = voe Ve 
without this term. In both cases the difference in the line posi- 
tion is a factor of e~ 2y ", a property that is in agreement with the 
more simple solution 

As an example, in Fig. [1] and [2] we show the time-evolution 
of an initially narrow line for different values of the y-parameter. 
For the curves shown in Fig. [T]the effect of stimulated scatter- 
ing in the blackbody ambient radiation field was included, while 
in Fig [2] it was not. As expected, in the former case the max- 
imum of the distribution is moving toward lower frequencies, 
while in the latter it moves toward higher frequencies. It is also 
remarkable that the analytic solution (fT4T > in comparison with 
the full numerical solution of Eq. (TT2l ) works up to values of 
y e ~ 1 - 5. It only breaks down when the photons reach the 
region x e > 1, where the recoil effect starts to become impor- 
tant. However, in both considered cases (injection at x e $ = 1CT 2 
and x e o = 0.1) only a small fraction of the initial number of 
photons reaches this region for y e < 1. One can also see that 
no matter if the background-induced effect is included or not, 
the line-broadening is not changing much. It is only the posi- 
tion of the line center that is affected by this process, but not the 
Doppler-broadening . 

5. Discussion and conclusion 

We have presented a new solution of the Kompaneets-equations 
for physical situations in which low frequency photons are scat- 
tered by thermal electrons but within an intense ambient black- 
body field that is very close to full thermodynamic equilib- 
rium with the free electrons. Under these circumstances the 
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Fig. 1. Time-evolution of AA^ V = v 2 An for different values of 
the y-parameter including the effect of stimulated scattering in 
the blackbody ambient radiation field. The upper panel shows 
the case, for an initially narrow line which was injected at fre- 
quency x e o = 10 2 , while the lower panel shows the solution 
for injection at x e o = 10 _1 . In both figures we show the results 
as obtained by numerically solving Eq. (TT~2T > for T y = T e . In ad- 
dition we give the analytic solutions of the linearized problem, 
Eq. (13[ , according to Eq. ( TBI ). 




Fig. 2. Time-evolution of AN V = v An for different values of 
the y-parameter not including the effect of stimulated scatter- 
ing in the blackbody ambient radiation field (n p i = 0). The 
upper panel shows the case, for an initially narrow line which 
was injected at x e o = 10~ 2 , while the lower panel shows the 
solution for injection at x e o = 10 _1 . In both figures we show 
the results as obtained by numerically solving Eq. (fTZb for 
T y = r e . In addition we give the analytic solutions according 
to lZeldovich & Sunvaevl fl969), Eq. ©. 



blackbody-induced stimulated scattering of photons slows down 
the motion of photons toward higher frequencies by a factor of 
3 in comparison with the solutio n that only take into accoun t 
Doppler-broadening and boosting (Zel dovich & Sunya ev 1969). 

These physical conditions are met during the epoch of 
cosmological recombination, in particular for He in — » He n- 
recombination, where electron scattering can still af- 
fect the cosmologica l reco m bination spectrum notabl y 
dDubrovich & Stolvarovl Il997t iRubino-Martin et alj |2007|) . 
The CMB blackbody spectrum is an exact solution of the 
Kompaneets equation when T e = T y . Nevertheless, the weak 
low frequency features in the presence of CMB with T e = T 7 
are evolving and boosting their central frequencies due to the 
simultaneous action of induced scattering and Doppler-effect 
resulting after each scattering. The additional line shift towards 
low frequencies due to recoil effect is small because x e <K 1 . 



Since at z ~ 6000, i.e. where most of the He n-photons are 
released, one finds y e ~ 2 x 10~ 3 , so that within the standard 
recombination epoch this process is only important for very 
accurate predictions of the cosmological recombination spec- 
trum at low frequencies. These may become necessary for the 
purpose of using the cosmological recombination spectrum for 
measurements of cosmological parameters, such as the CMB 
temperature monopole, the specific entropy of the Universe, 
or the pre-stellar abunda nce of helium, as discussed recently 
dSunvaev & Chlubal 120081) . Our analysis shows that the low- 
frequency spectral features from He in — > He II recombination 
should then be shifted by Av/v ~ 0.2% instead of ~ 0.6%. 
At the same time the Doppler-broadening of these lines reaches 
-15% at FWHM, and the t ypical line width is Av/v ~ 20 - 30% 
(IRubino-Martin et al.ll2007l) . 

However, if photons due to uncompensated atomic transi- 
tions in hydrogen and helium were released much earlier (z > 
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10 4 - 10 5 ) then the blackbody-induced stimulated scattering be- 
comes more significant. As mentioned in the introduction, such 
emission can occur in connection with possible intrinsic spec- 
tral distor tions of the CMB, due t o energy release in the early 
Universe dChluba & Sunvaevl2008l) . 

We should stress that the discussed behavior of spectral fea- 
tures is characteristic only for the low frequency part of CMB 
spectrum hv <K kT y . In the dm and cm spectral range this con- 
dition is fulfilled with very good accuracy, and the occupation 
number of photons exceeds ~ 10. The broadening and shifting 
of the photons in the Wien-part of the spectrum has completely 
different behavior and will be discussed in a separate paper. 

In addition we gave the solution (|6]l which describes the evo- 
lution of narrow spectral features in the presence of a very bright, 
broad-band radiation background with «back » 1 . It might have 
interesting applications for radio lines arising inside of extra- 
galactic radio source, where according to our estimates (see 
Sect. [3]) the background-induced line shifts could be amplified 
and even exceed the line broadening due to the Doppler-effect. 
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